Solid Surface Structure Affects Liquid Order at the Polystyrene/SAM Interface 
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We present a combined x-ray and neutron reflectivity study characterizing the interface between 
polystyrene (PS) and silanized surfaces. Motivated by the large difference in slip velocity of PS on top 
of dodecyl-trichlorosilane (DTS) and octadecyl-trichlorosilane (OTS) found in previous studies, these 
two systems were chosen for the present investigation. The results reveal the molecular conformation 
of PS on silanized silicon. Differences in the molecular tilt of OTS and DTS are replicated by the 
adjacent phenyl rings of the PS. We discuss our findings in terms of a potential link between the 
microscopic interfacial structure and dynamic properties of polymeric liquids at interfaces. 



I. INTRODUCTION 

When downsizing devices, confinement and interface 
effects grow enormously in importance. Apart from the 
fundamental interest in interfacial structure and dynam- 
ics, micro- and nanodevice fabrication open new perspec- 
tives for applications in manufacturing, pharmaceutics, 
chemistry or the food industry [TH3]. Especially, in the 
context of microfluidic devices £Q, the controlled mo- 
tion of small amounts of liquid is indispensable. As the 
solid/liquid friction dramatically impacts hydrodynamics 
in these systems, the boundary condition (BC) of flowing 
liquids, commonly quantified by the slip length 0], has 
been extensively revised on a microscopic length scale in 
recent years [5]. 

Navier first supposed that a liquid may slip over a solid 
surface [4] and introduced the slip length 6, which is de- 
fined as the distance z from the interface where the ve- 
locity profile v(z) of the liquid extrapolates to zero: 

b = v(0)(dv(z)/dz)- 1 \^ . (1) 

The possibility of surface slippage was subsequently 
intensely discussed. Since the mid 19th century, the 
no-slip BC had been generally favored and in the 20th 
century, fluid dynamics textbooks assumed it in general, 
often without reference to its empirical origin. Experi- 
mental evidence of the failure of the no-slip BC was often 
attributed to parasitic effects or the lack of resolution, 
until de Gennes theoretically predicted large slip of 
entangled polymers in capillaries [6]. This provoked 
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many investigations and today the slippage of entangled 
polymer melts is a well-known phenomenon [7]. In the 
late 1990s, slip of Newtonian liquids was observed in 
molecular dynamics (MD) simulations [5] and shortly 
after that it was experimentally confirmed [HI fTO] . 
Although interfacial slip developed to a well-recognized 
phenomenon, its microscopic origin is still unclear, also 
hindered by the fact that most experimental techniques 
used to determine the slip length are invasive or indirect. 
From a theoretical point of view, two different types of 
slippage are distinguished [11]: Real slip occurs when 
the liquid slides over the solid surface on an atomic scale. 
Alternatively, apparent slip arises where a microscopic 
boundary layer exists that is structurally and/or dy- 
namically different from the bulk liquid. This boundary 
layer may lead to a different viscosity and is observed 
as interfacial slip on a larger length scale, although the 
no-slip BC may microscopically still hold. The nature of 
such a boundary layer may be a depleted density of the 
liquid [12 or an alignment of the near-surface molecules 

Depletion effects of simple liquids have been observed in 
various cases |14H17j using x-ray and neutron reflectom- 
etry (XRR and NR). Their origin and the consequential 
link to macroscopic properties of liquids on solid surfaces 
such as hydrophobicity and also slippage is currently 
under debate [TtjHTB] . Density profile fluctuations 
and depletion layers of polymer melts close to solid 
substrates were also reported and attributed to altered 
molecular conformations and locally modified segmental 
distributions [19] . In cases of entangled polymer melts, 
dedicated chain conformations at the solid/liquid inter- 
face are responsible for a decrease in the entanglement 
density compared to the bulk and, thus, substantially 
influence slippage [217] . 

In this work, we present a combined XRR and NR 
study on polystyrene (PS) films on top of two dif- 
ferent silanes, octadecyl-trichlorosilane (OTS) and 
dodecyl-trichlorosilane (DTS). As known from previous 
studies, PS melts show large slippage when flowing 
over hydrophobized surfaces, depending on molecular 
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weight [2U], temperature and substrate [2T]. In contact 
with DTS, a slip length of roughly 1 fj,m (cf. Ref. [2"T] ) 
has been found for PS of the same molecular weight 
(13.7kg/mol) and for the same annealing temperature 
(120 °C) used in this study. On OTS, however, the 
sip length is about one order of magnitude shorter for 
the same parameters. This strong effect on the slip 
length is surprising, as both are chemically identical 
self-assembled monolayers (SAM) that differ only by 
six backbone hydrocarbons in tail length. We show 
that the difference between the two surfaces is a lower 
grafting density of the DTS resulting in a tilt of the 
hydrocarbon tails. We provide evidence that the SAMs 
induce conformational changes within the interfacial 
polymer, which may influence slippage. 



II. EXPERIMENTAL SECTION 

The Si(100) wafers (Wacker/Siltronic, Burghausen, 
Germany, boron p-dotation, 10-20 f2cm resistance) were 
hydrophobized with OTS and DTS monolayers [22] . 
which resulted in a static contact angle of 67 ± 3° for 
PS on the silanized wafers in both cases. The advancing 
water (Milli-Q synthesis system, Millipore, USA, organic 
impurities <6ppb, resistance at 25°C: 18.2 MSlcm) 
contact angle was 116° on OTS and 114° on DTS. 
The receding contact angle was 110° in both cases. 
The atactic PS with a molecular weight of 13.7kg/mol 
(M w /M n =1.03) and the deuterated PS (dPS) with a 
molecular weight of 12.3kg/mol (M w /M n =l.Q5) were 
purchased from PSS, Mainz, Germany. PS films between 
50 nm and 60 nm were prepared by spin-casting a toluene 
(Merck, Darmstadt, Germany) solution onto mica and 
floating on Millipore water, from where they were picked 
up by the hydrophobized wafers. Then, the samples 
dedicated for the reflectivity experiments were annealed 
above the glass transition temperature (T g ) at 120 °C 
for 30 s just before the onset of dewetting. Results of the 
dewetting studies and the slip length determination can 
be found elsewhere [21j . The bare silanized substrates 
were measured unannealed. 

The x-ray measurements were conducted at beam- 
line BL9 [23 of the Dortmund Electron Accelerator 
(DELTA), Germany, with photon energies of llkeV 
and 15.2 keV and beam sizes of 0.2*2.5 mm 2 and 
0.1*1 mm 2 , respectively, with an angular resolution 
of 0.008° (FWHM). We observed no beam damage 
during the x-ray measurements. The NR measurements 
were performed on the ADAM Reflectometer [21] at 
the Institut Laue-Langevin (ILL) in Grenoble, France, 
using a 0.5*10 mm 2 beam with a constant angular 
resolution of 0.08° (FWHM). Pre-characterization of the 
samples dedicated for neutron scattering by means of 
XRR was done on the laboratory reflectometer XPcrt 
Pro PW3020, Panalytical, Netherlands. Fitting of the 
reflectivity data was obtained by using co-refinement of 



a slab model with Motofit [25] . The error bars given by 
the fit correspond to one standard deviation, whereas 
the parameter ranges indicated in this work correspond 
to the variation among two measurements of equal 
samples from different batches. 



III. RESULTS 
A. X-ray reflectivity of bare silane substrates 

The x-ray reflectivities normalized to the Fresnel 
curve, i.e. the reflectivity for an ideally flat silicon 
surface, are shown in Fig. 



1(a) and (b). To get a 



quantitative description of the measured results, we first 
analyze the silanized Si wafers and assume a three-slab 
model, consisting of silicon oxide (Si02), a silane 
headgroup and a hydrocarbon tail [25]. This reproduces 
the measured data, as can be seen by the solid lines 
in Fig. [I] The SAM parameters (see supplementary 
information) reflect the characteristics of completely 
grown silane layers [25] [27]: On top of the 9- 10 A thick 
Si02, with densities of 2.24- 2.25 g/cm 3 , there is a silane 
headgroup, which is 5.6 -5.95 A thick, followed by the 
silane tail. The roughness between subsequent layers 
is between 1-4 A. Note that the roughness of the bare 
DTS tail of 2.9 ± 0.08 A is comparable to the roughness 
of 2.73 ± 0.01 A measured for OTS. This rules out the 
possibility of surface roughness as a potential cause for 
difference in slippage between OTS and DTS. Moreover, 
in-plane rocking curves at different q z values on the 
reflectivity showed no significant broadening of the spec- 
ular peak on both silanized surfaces which underlines 
the flatness and homogeneity of both samples. The only 
significant difference between the DTS and the OTS 
layer, apart from the tail length, is a higher grafting 
density of the OTS, which appears as a higher electron 
density of the silane head and tail. The headgroup of 
the DTS has an electron density of 0.476 ± 0.007 A~ 3 , 
compared to 0.532 ± 0.004 A" 3 in case of the OTS. 
Likewise, the density of the DTS tail (0.82 ± 0.01 g/cm 3 ) 
corresponds to 88% of an alkane crystal's density [28], 
whereas the OTS reaches 100% (0.936 ± 0.004g/cm 3 ). 
This difference is also observed when comparing the 
thicknesses of the layers with the calculated length of 
an all-trans hydrocarbon chain [25] . The 21.31 ± 0.05 A 
OTS tail length matches 99 % of the calculated fully 
stretched molecule (21.5 A). The 12.0 ± 0.1 A DTS 
tail length, however, corresponds to only 86 % of the 
calculated 13.9 A, which is commonly explained as a 
randomly tilted SAM [25J [23] and yields a tilt angle 
of 30° for DTS in this case. This difference between 
DTS and OTS is regularly observed [30] which may be 
due to the less optimal preparation temperature (room 
temperature) of DTS in comparison to the OTS [21] ■ 
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Figure 1 . X-ray reflectivity curves normalized to the Fresnel-reflectivity on a logarithmic scale of the bare silicon wafers covered 
with DTS (purple bars) (a) and OTS (green bars) (b) as well as the reflectivities for the silanized samples covered with PS 
(shifted by one, red and blue bars, respectively). The solid lines represent fits. Fig. (c) and (d) display the corresponding 
scattering length density (SLD) profiles in the same color code as the reflectivity data. To exclude interfacial roughness as the 
origin of the depletion layer, the dotted black line (d) was obtained from adding the SLDs of the OTS and PS profiles in air. 
The resulting simulated reflectivity in (b) demonstrates that the experimental reflectivity data cannot be explained by just a 
roughness mismatch between OTS and PS. In order to highlight the difference in the depletion layer thickness for OTS and 
DTS, the blue dash-dotted lines in (a) and (c) denote a simulation of the DTS data with a thicker interfacial layer. 



B. X-ray reflectivity of PS on silane substrates 



When the silanized substrates are brought into contact 
with PS and annealed well above the bulk T g of PS, the 
data analysis reveals an interface layer with lower density 
between the SAM tail and the PS as can be seen in Fig. 
[l] To check whether this density-reduced layer is just an 
artifact of a roughness mismatch of the adjacent layers 
due to, for example, insufficient annealing we calculated 
the electron density profile of the interface between 
the OTS and PS layers by just adding both profiles 
measured in air before contact, assuming no change in 
the layers themselves. The simulated scattering length 
density (SLD) profile and the resulting reflectivity are 
denoted by black dotted lines in Fig. [ijb) and (d). This 
assumption does not reproduce the measured curve, 
obviously the PS has changed in contact with the SAM. 
The SLD profile corresponding to the best fit to the 



reflectivity data which is shown by the solid lines in 
the same figure, features a sharp density change of the 
interface layer towards the residual PS film which points 
out that a smooth interface has emerged on a molecular 
level. The relevant fit parameters are summarized in 
Tab. [T] in the supplementary information. The major 
difference between the DTS and OTS interfaces is the 
thickness of the density-reduced interface layer. The 
4.2 ± 0.14A thick layer at the PS/DTS interface is 
considerably thinner than the 6.79 ± 0.04 A at the 
PS/OTS interface. 



C. Neutron reflectivity of PS on silane substrates 

In order to obtain information about the chemical 
composition of the interface layer between OTS/DTS 
and PS, we have performed NR experiments on the 
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Figure 2. Neutron reflectivities in log scale of the OTS (blue 
bars, shifted by one) and DTS (red bars) substrates covered 
with deuterated PS. The solid lines are fits corresponding to 
the SLD profile from the inset in the same color code as the 
data points. 



same system. In contrast to x-rays, which are sensitive 
to the electron density of the sample, neutrons are 
scattered by nuclei and the scattering length difference 
between a proton (-3.7 fm) and a deuteron (6.7 fm) 
is noticable. This makes NR particularly sensitive to 
protonated / deuterated interfaces. 

Replacing the PS by dPS, we obtained the NRs shown 
in Fig. [2j As the neutron measurements suffer from 
a considerably smaller momentum transfer range as 
compared to XRR their spatial resolution is much lower. 
Thus we fixed all fit parameters when analyzing the NR 
measurements by the corresponding x-ray values of the 
hydrogenated samples and varied only the SLD of the 
interfacial layer (and the SLD and thickness of the PS 
layer). For those parameters neutrons are more sensitive 
than x-rays due to the large contrast between protonated 
silane (SLD = -0.4 * 1(T 6 A- 2 ) and deuterated PS (SLD 
= 6.6*KT 6 A- 2 ) (for x-rays the contrast is about 10 
times smaller). This gives us a clear picture of the 
isotope composition of the interfacial layer. A purely 
protonated silane layer would result in a small negative 
SLD whereas a deuterated PS layer would show up in a 
layer with a SLD of about 4.5 * lO^A" 2 . 
In the inset of Fig. [2] it is evident that none of the 
two scenarios is observed for the SLD of the interfacial 
layer. Instead both interfaces seem to consist of both 
protonated and deuterated material. By taking the 
electron density from the x-ray SLD and the nuclear 
density of the NR SLD the exact amount of PS and 
Silane of the density-reduced layer can be calculated. 
For the OTS interface, this results in 66 ± 4 % PS and 
11 ± 3 % silane as compared to their bulk density, and, in 
case of the DTS, 32 ± 12 % PS and 43 ± 12 % silane are 
present. This means that the observed low-density layer 
in the x-ray measurements comprises parts of the SAM 
and the adjacent PS. This result is in-line with recent 
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Figure 3. Sketch of the molecular conformation at the 
PS/OTS (a) and PS/DTS (b) interface as explained in the 
text. 



x-ray reflectometry studies of water at hydrophobic 
surfaces [T5J [TB] and clarifies that the density-depleted 
liquid often observed at solid/liquid interfaces is partly 
due to the presence of the hydrogen termination of the 
hydrophobic SAM whose SLD is close to for x-rays. 



IV. DISCUSSION 

To produce the sharp density step of PS in contact 
with the SAMs, as seen by the electron density profile, 
the adjacent PS chains cannot be randomly oriented. In- 
stead, we assume that a rather well-ordered arrangement 
of contacting chain segments is formed. Considering the 
molecular composition of PS, only an orientation with 
the phenyl groups pointing to the SAM complies with all 
parameters extracted from the scattering experiments. 

The distance 



This scenario is sketched in Fig 3(a) 
between the hydrocarbon backbone of the PS and the 
end of the phenyl group, including the covalent radius 
of the hydrogen, is 5.6 A. The projected bonding length 
of the OTS hydrogen termination including its covalent 
radius is 0.65- 1.25 A, depending on whether the covalent 
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radius of the linked carbon is subtracted or not. In total 
this adds up to 6.25 - 6.85 A, in accordance with the 6.8 A 
thick interface layer observed at the PS/OTS interface. 
Due to three possible orientations of the phenyl group 
around a flat PS backbone, in average only every third 
one would be incorporated in the interface layer. In 
the bulk, one and a half out of three phenyl rings are 
projected on one side of the PS backbone. This explains 
the density of deuterated material, reduced by one third, 
in the interface layer revealed by NR. If we additionally 
assume that the phenyl rings follow the orientation of 
the SAM , they would be tilted by 30° as depicted in 
Fig. |3(b)| in contact with DTS. This would lead to a 
reduced interface layer of 4.18- 4.78 A, which matches 
very well with the 4.2 A thick interface layer observed 
at DTS. In contrast to one methyl- hydrogen pointing 
into the interface layer at OTS, two methyl- hydrogens 
are present in the interface layer of the tilted DTS. 
Additionally, the interface layer at DTS is thinner than 
the OTS one and hence, the proportional amount of 
protonated silane should be considerably higher at the 
DTS, at the expense of the relative amount of deuterated 
PS. This is well in accord with the 43 % silane and 32 % 
dPS as revealed by NR for the DTS interface. 
The interfacial order which we deduce from our scatter- 
ing experiments is in line with recent MD simulations, 
where a crystalline surface was able to induce order in 
polymeric liquids |32j : The first liquid layer showed 
an almost perfect reproduction of the preset periodic 
crystal structure. Recent experiments confirm, that in 
thin polymer films chain segments may order [33] and, 
in particular, that certain orientations of the PS phenyl 
rings can be induced by the presence of interfaces [34], 
even in absence of specific interactions: As demonstrated 
by non-linear optical techniques such as sum-frequency 
generation (SFG) spectroscopy, the interplay of intra- 
and intermolecular interactions causes the phenyl rings 
to point away from the bulk polymer film perpendicular 
to the polymer/air interface |35j and also towards a 
hydrophobic substrate |36j . 

Concerning slippage, MD studies report slip lengths 
on the order of several monomer lengths [32 which 
are quantitatively not comparable to the large experi- 
mental slip lengths (up to several fim close to T g ) for 
unentangled polymer melts on silanized surfaces |21j . 
Experimentally determined slip lengths may incorporate 
apparent slip in addition to real slip. Apparent slip 
of polymeric liquids may be explained by a higher 
segmental mobility in the vicinity of the interface, either 
(i) due to an interfacial depletion effect [T3] or (ii) due 
to layering [SJ |SZ] and/or alignment [T3] of the liquid 
near the interface. As a third potential mechanism 
(iii), a reduced segmental friction coefficient between 
the adjacent polymer chains and the substrate may 
be caused by particular polymer conformations near 
interfaces. Recent MD simulations on slipping oligomers 
[55] highlight that the most significant part of the energy 
transfer (friction) between the solid and the liquid is 



dissipated in the first liquid layer and only minor energy 
transfers occur between subsequent layers. 
The structural data of our study imply a flat arrange- 
ment of the adjacent PS chain segments and a sharp 
step in the density profile between the interfacial layer 
and the residual polymer film. The density was shown 
to be reduced down to 75 - 77% in a depletion zone 
of 4 - 7 A thickness. Mechanism (i) does not apply to 
our system since the small extent of the depletion zone 
cannot account for the large experimental slip lengths 
measured. Regarding mechanism (ii), apparent slip due 
to layering implies more than one layer of polymer to be 
aligned. However, density oscillations indicating molec- 
ular layering were not detected. Such shear-induced 
effects have been observed in in situ experiments for 
a different system while applying shear flow 39J. Our 
experiments clearly demonstrate a distinct orientation 
of the phenyl rings of the PS melt due to the structural 
properties of the adjacent substrate (OTS and DTS). 
Hence, compared to the non-oriented bulk liquid, the 
presence of locally deviating dynamical properties such 
as friction and viscosity in the interfacial region are very 
likely. 

Theoretical studies report a suppression of slip for a 
higher degree of substrate-induced liquid ordering near 
a flat smooth surface (c.f. [1Q] and referenes therein). 
This is qualitatively in line with the difference in slip 
length that has been experimentally observed for PS 
on OTS and DTS. In general, an ultra-low surface 
roughness of the substrate, as it is found in the case 
of OTS and DTS, is a precondition for the presence 
of ordering effects. In the limit of commensurability 
of spatial dimensions of the surface corrugations and 
the molecular size of the liquid, a strong suppression of 
slip can be achieved [10]. This might explain small slip 
length values (between and 100 nm at maximum) that 
have been experimentally found for the same PS on an 
amorphous substrate exhibiting slightly larger surface 
roughness [2U] . Studies of the interface structure of the 
latter system are on-going. 

In literature, experimentally observed differences in slip 
length of Newtonian liquids on different substrates have 
been widely attributed to surface properties such as 
roughness and the strength of interaction between liquid 
molecules and the substrate [HI SI] ■ We stress the fact 
that these parameters are found to be identical for PS 
on DTS and on OTS. For these systems, we provide 
evidence of a molecular interplay of the interfacial 
structure of the liquid and the surface order of the 
solid, that might affect a macroscopically detectable 
parameter, namely the amount of slip of a PS film on 
silanized surfaces. 
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V. SUMMARY AND OUTLOOK 

In summary, we have revealed that surface structure 
of a self-assembled monolayer affects the conformation of 
polymer chain segments adjacent to the solid boundary. 
The results of combined x-ray and neutron scattering 
studies point out that the adjacent polystyrene chain 
segments lie completely flat and, moreover, that the ori- 
entation of the phenyl rings replicates the self-assembled 
monolayer structure. Both facts appear to be the 
clue to understand a) substantially different polymer 
slippage on silanized surfaces exhibiting identical surface 
energies and polystyrene contact angles and b) large 
effective (comprising real and apparent) slip, both 
observed experimentally. Our findings corroborate 
on-going research claiming conformational changes at 
the interface in case of entangled polymeric liquids 
|20j and the interfacial liquid structure in case of non- 
cntangled oligomers as important parameters governing 
macroscopic slip [T71 02] • Additionally, our results might 



also shed light on further interfacial phenomena such 
as depletion layers or glass-transition temperatures of 
thin polymer films. Since the orientation of the phenyl 
rings of the polystyrene is linked to the aforementioned 
ordering phenomena, molecular dynamics simulations 
investigating slippage of polymer melts should intend to 
account for their entire monomeric structure to achieve 
full comparability to experimental situations. 
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SUPPLEMENTARY 



Table I. Fitting parameters of the silane tail and the interfacial layer at the SAM/PS interface measured with X-ray reflectometry. 
d, Nb and a correspond to the thickness, scattering length density and roughness of the respective layers. 



Parameter 


DTS in air 


OTS in air 


PS on DTS PS on OTS 


ivf [nr 6 A- 2 ] 


20 ±0.1 


20 ±0.1 


20 ±0.1 


20 ±0.1 


aSl /sw 2 [A] 


3.3 ±0.8 


1.1 ±0.1 


3.3 ±0.8 


1.1 ±0.1 


d Si °* [A] 


9± 1 


10 ±0.1 


9± 1 


10 ±0.1 


N Si0 2 [ 10 -6^-2j 


18.8 ±0.3 


19.08 ±0.05 


18.8 ±0.3 


19.08 ±0.05 


o Si0 2 /head J^J 


2.1 ±0.5 


1.4 ±0.1 


2.1 ±0.5 


1.4 ±0.1 


d head ^ 


5.93 ±0.02 


5.62 ±0.02 


5.93 ±0.02 


5.62 ±0.02 


N head [1Q-6A- 2 ] 


13.4 ±0.2 


15 ±0.1 


13.4 ±0.2 


15 ±0.1 


^.head/tail rj^j 


2.5 ±0.2 


2.19 ±0.04 


2.5 ±0.2 


2.19 ±0.04 


d tail 


12 ±0.1 


21.31 ± 0.05 


8.84 ±0.08 


17.86 ±0.03 


Nl ai1 [io~ 6 A- 2 ] 


7.9 ±0.1 


9.07 ±0.02 


7.9 ±0.1 


9.07 ±0.02 


ail 1 air j^j 


2.9 ±0.08 


2.73 ±0.01 






^t ail / inter f jj^j 






1.14 (fixed) 


1.14 ±0.03 


d interf [A] 






4.19 ±0.14 


6.79 ±0.04 








5.99 ±0.03 


6.44 ±0.03 


G interf/PS J^J 






(fixed) 


0±0.2 


d PS [A] 






539.9 ±0.6 


615.3 ±0.5 


N b PS [io- 6 A- 2 ] 






8.2 ±0.5 


8.2 ±0.5 


a PS/air [A ] 






3.2 ±0.02 


3.31 ±0.04 



